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Abstract. Tsallis Statistics was used to investigate the non-Boltzmann distribution 
of particle spectra and their dependence on particle species and beam energy in the 
relativistic heavy-ion collisions at SPS and RHIC. Produced particles are assumed 
to acquire radial flow and be of non-extensive statistics at freeze-out. J/-0 and the 
particles containing strangeness were examined separately to study their radial flow 
and freeze-out. We found that the strange hadrons approach equilibrium quickly from 
peripheral to central A+A collisions and they tend to decouple earlier from the system 
than the light hadrons but with the same final radial flow. These results provide 
an alternative picture of freeze-outs: a thermalized system is produced at partonic 
phase; the hadronic scattering at later stage is not enough to maintain the system 
in equilibrium and does not increase the radial flow of the copiously produced light 
hadrons. The J/-0 in Pb+Pb collisions at SPS is consistent with early decoupling and 
obtains little radial flow. The J/ip spectra at RHIC are also inconsistent with the bulk 
flow profile. 
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1. Introduction 

Identified particle spectra in transverse momenta are one of tlie pillars in the major 
discoveries in high-energy nuclear physics [11 [2l [3l H]. It has been demonstrated that 
the spectral shape is sensitive to the dynamics of the nucleus-nucleus collisions [5l E] 
and can be used to obtain the radial flow and temperature at freeze-out. In addition, 
it has been argued that hadrons containing strange or charm valence quarks should 
have smaller hadronic interaction cross-section, and should decouple from the system 
earlier than the hadrons with only light valence quarks [11 |5l [6]. In this scenario, those 
strange and charmed hadrons would carry direct information from the collisions at earlier 
stage without dilution due to hadronic scattering at late stage. Particle spectra at SPS 
and RHIC have been used to extract this behavior [5l [6]. The Blast-wave model with 
Boltzmann-Gibbs (BGBW) distribution has been applied to SPS and RHIC, and was 
the basis for the observation of the early decoupling of multi-strange hadrons [3 16]. On 
the other hand, the evolution with hydrodynamics shows that the multistrange particle 
spectra can be well described by the same hydrodynamics at the same freeze-out as 
other hadrons [7|. 

Recently, non-extensive hydrodynamics has been proposed to explain some essential 
features in relativistic heavy-ion coUisions [8l|9l[l0l[IIl[T2l[T3l[Tl[l5l[T6l[n]. A simplified 
model with blast-wave assumption and non-extensive Tsallis statistics for hadrons at 
freeze-out has been implemented and applied successfully to describe the RHIC data 
including all the multistrange hadrons upto intermediate [ID]- Same approach was 
extended to fit the hadron spectra at SPS [18]. Although Ref. [10] concluded that 
a common non-extensive Tsallis Blast- Wave model (TBW) can describe hadrons with 
and without strange valence quarks, Ref. [18] claimed that multistrange hadrons and 
^/ip freeze out at much higher temperature and acquire smaller radial flow than light 
hadrons. 

In this paper, we perform a systematic study of hadron spectra at SPS and RHIC 
top energies with TBW as in [10]. Detailed comparisons of spectral shape from model 
and data are carried out to assess the degree of this discrepancy. In addition, we 
categorize the spectra into all hadrons, strange hadrons, non-strange hadrons and ^/ip 
in an attempt to compare their spectral shapes to the TBW results. 

2. TBW and its fit to RHIC and SPS data 

The same TBW as in Ref. [10] is used in the current analyses: 



where the left-hand side is invariant differential particle yield at mid-rapidity, mx 




Pt sinh(p) cos(0))) 



(1) 
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and Pt are transverse mass and transverse momentum of the produced particle, q is the 
parameter characterizing the degree of non-equihbrium, and p is the flow profile growing 
as n-th power (n = 1) along the transverse radial direction (r) from zero at the center 
of the collisions to /3s at the hard-spherical edge (i?). 

The STAR Collaboration has published a series of particle spectra at mid-rapidity. 
The most complete set is for p-|-p and Au+Au collisions at y/s^ = 200 GeV. The 
identified particle spectra include tt^, K"^, Ks, K*, p, (p, A, H, Q, p, A, H, and 
0[Il[I21l2niEIll221l231[2ai2ai2Sl[27]. In addition, the PHENIX Collaboration has 
published spectra of charged kaons [28] and J/xp [291 EOl EI], which are not included in 
the previous fits [TD] . The charged kaon spectra from PHENIX Collaboration extend the 
STAR measurement to higher p^. This is crucial for the current study when the spectra 
are categorized into several groups. In the group of hadrons with strange valence quarks, 
kaon is the only light meson, whose spectrum is more sensitive to the non-equilibrium 
effect at low since the spectra of heavier particles are usually overwhelmed by radial 
flow in the low-px to intermediate p^ in the central nucleus-nucleus collisions. J/ip 
spectra are measured by PHENIX Collaboration inp+p, Cu+Cu and Au+Au collisions 
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Figure 1. (Color Online) Identified particle spectra from Au+Au collisions at ^s„„ 
=200 GeV. The left panel: spectra of light hadrons and strange hadrons. The solid 
curves are results from TBW fit. the right panel: J/V' spectrum, the dashed line is the 
TBW prediction using parameters from the fit to other hadrons, and the solid curve 
is a TBW fit to 3/ip alone 
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Figure 2. (Color Online) Identified particle spectra in central Pb+Pb collisions at 
beam energy of 158 AGeV from the fixed-target experiments at SPS. The left panel: 
spectra of light hadrons and strange hadrons. The solid curves are results from TBW 
fit. The right panel: J/ip spectra, the dashed line is the TBW prediction using 
parameters from the fit to other hadrons, and the solid curve is a TBW fit to 3/4' 
alone. 



at RHIC top energy for pt < 5 GeV/c [29l EQl EI] while STAR Collaboration has 
extended these measurements to high-pT [32] • Only PHENIX data with < 5 GeV/c 
are used in this TBW analysis. We perform a TBW fit using ([1]). The procedure and 
program are described in details in [101 [33] . 

Figure [1] shows the identified spectra measured by STAR and PHENIX 
collaborations and their associated TBW curves from the fit. The yields are presented 
in terms of invariant differential cross-section. There are a few features: the extension 
of kaon spectra to high p^ doesn't change the quality of the fit and the kaon spectra 
can be well described by the TBW. However, using the same parameters from the fit, 
the TBW curve can not reproduce the J/ifj spectra as shown in the right panel of 
figure [H The curve is too fiat at low p^, presumably due to the effect of large radial 
flow. We performed a TBW fit to the J ftp spectra alone. The result shows that J ftp 
has consistently lower radial flow than the bulk, although the parameter uncertainty is 
quite large due to the large statistical uncertainties of the data points. 

Similarly, the same TBW fit was applied to SPS data [3ll[35l[36l[371[38l[39l[lQl[lIl 




Figure 3. (Color Online) 1-a contour of T vs (q— 1) plane from TBW fit to the J/ip 
spectrum in p+p collisions at y/s =200 GeV with x^/nDoF=13/14 at minimum (a), 
Pb+Pb collisions at =17.2 GeV with xVnDoF=153/38 (b) and ratio of data 

points and fit in Pb+Pb collisions (c). The radial flow is set to zero in this fit. 



W2\ . Figure H] shows the measured spectra together with the fit TBW curves for central 
Pb+Pb and In+In collisions at ^/s^ =17.2 GeV. The TBW results are generally in 
good agreement with the data points. However, the x^/nDoF is not as good as the 
fit at RHIC. Although the TBW curves represent the trends of the data points very 
well, the variations of the data points around the curves seem to be larger than those 
allowed by the error bars, indicating a possible underestimate of the uncertainties. In 
right panel of figure [2l the J/ijj data points from NA50 and NA60 were presented as a 
function of together with the TBW curves. The dashed line is the prediction using 
the parameters produced by the TBW fit to the other light hadrons, while the solid 
curve is a TBW fit to the J/ip spectrum alone. Same situation as observed at RHIC, 
the TBW prediction with parameters obtained from the fit to the light hadron spectra 
overpredicts the radial flow present in the J/ip spectrum. A fit to J/?/' spectrum alone 
shows that the radial flow of J/ip at SPS is consistent with zero. 

To further investigate the effect of the correlation between T and g as a potential 
cause for the large uncertainties from the J/ip fit, we examine the contour plot from the 
TBW fit to the J/ip spectra. The J ftp spectrum in p+p collisions at ^/s =200 GeV has 
much smaller uncertainties with many more data points than those in A+A collisions 
at RHIC. We perform a fit by setting radial flow to zero while only keeping T and q 
as free parameters. Figure |3] shows the 1-a contour of the T vs (g — 1) plane with the 
best fit as presented in figure [1] and figure [2l It is obvious from the plot that there is 
a strong correlation between these two parameters. Once one parameter is constrained 
accurately by other means or other hadron spectra, the other parameter is also well 
constrained. 

Intuitively, one expects that the strange hadrons with smaller hadronic interaction 
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Figure 4. (Color Online) q, T, /3 and x^/nDoF as a function of centrality for different 
groups of hadrons from the TBW model fit to spectra in Au+Au collisions at RHIC 
top energy. 



Table 1. Fit results for all hadrons, strange hadrons and non-strange hadrons in 
central A+A colhsions at RHIC and SPS. 
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cross-sections would freeze-out earlier than the bulk of hadrons [5] . This results in higher 
temperature and smaller radial flow for strange hadrons than for light hadrons without 
strange content. To study this effect, the particles are further grouped according to 
their species: all hadrons, strange hadrons and non-strange hadrons. Table [1] shows 
the fit results for these groups for central nucleus-nucleus collisions at SPS and RHIC. 
Overall, reasonable x^/nDoF is present in all cases. We note that the fits to the strange 
hadrons produce small x^/nDoF (~0.5), indicating an experimental overestimate of the 
systematic uncertainty. The strange hadrons consistently show smaller non-equilibrium 
effect as manifested in a smaller (q — l) value with higher temperature and similar radial 
flow than the non-strange hadrons. On the other hand, the TBW fit to the strange 
hadrons at SPS shows larger value of x^/nDOF, indicating less than ideal description 
of the data. In addition, the Q and Cl data points have larger fluctuations around 
the smooth TBW lines than the error bars allow. Figure |4] presents the q, T and P 
dependence on the centrality in Au+Au collisions at RHIC top energy. In fact, the 
(g — 1) is consistent with zero for strange hadron group and is significantly smaller than 
that of non-strange hadrons in the 0-40% central Au+Au collisions at RHIC while they 
approach each other in peripheral bins (40-80%). This is consistent with the findings of 
Refs. [261 Il2] • the strangeness saturation factor 7., increases from peripheral to central 
A+A collisions towards unity at RHIC and SPS top energies, indicating strangeness 
equilibration in central A+A collisions at RHIC and SPS. 

3. Discussions 

It is clear that the TBW model can provide satisfactory description of the spectra 
in nucleus-nucleus collisions at RHIC and SPS energies. There are several important 
conclusions from the results: 

(i) The radial flow velocity at SPS is smaller than that at RHIC. 

(ii) Freeze-out temperatures for non-strange hadrons are similar at RHIC and SPS. 

(iii) The non-equilibrium parameter (q — l) is small in central nucleus- nucleus collisions 
at RHIC and SPS except a larger (g — 1) value for non-strange hadrons at RHIC 
energy. 

(iv) The categorized groups of strange and non-strange hadrons show different 
parameters, indicating a possible early freeze-out of strange hadrons at RHIC. 
The biggest difference for strange hadrons from the light hadrons is the freeze- 
out temperature: strange hadrons freeze-out at higher temperature with similar 
radial flow. A possible conclusion is that the hadronic phase doesn't increase the 
radial flow of light hadrons at RHIC energies. The picture at SPS is sort of unclear: 
Although the strange hadrons have slightly higher temperature, smaller radial flow, 
and smaller (q-1) than those of the non-strange hadrons, the difference is at about 
2a level. In addition, the pr reach for the strange hadrons is considerably lower 
than those for the non-strange hadrons and the fit quality is less than ideal. 
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(v) The centrality dependence of (g — 1) of strange hadrons shows that the strange 
hadrons are off-equihbrium in peripheral colhsions as other hadrons. However, 
they approach thermal equilibrium very quickly at middle centrality with very 
small (g — 1) value. 

(vi) Most importantly, the J/?/) spectra at SPS and RHIC require much smaller flow 
velocity, and much larger off-equilibrium and/or higher temperature than the 
other hadrons. It also shows that the spectral shapes are the same in p+p and 
central nucleus-nucleus collisions. This is inconsistent with thermal production 
of J ftp through recombination of thermal open charms at chemical freeze-out in 
the central Au+Au collisions at RHIC [HI HSl |l6l S^. On the other hand, a 
dynamic coalescence throughout the system evolution [IH] concludes that J/ip 
acquires little flow in the process. This model doesn't maintain J/ip in equilibrium 
and is consistent with our conclusion. 

We noted that Ref. [H] performed a combined fit to the strange hadrons and J/ip , 
and concluded that both multi-strange hadrons and J/ip freeze-out early. However, our 
separate fits provide a different picture, where multi-strange hadrons show deviations 
from the bulk hadrons with high value of temperature while J/ip spectral shape shows 
little bulk effect, and they are incompatible with each other. The forced combination of 
these two groups doesn't provide a good fit, but does produce a set of parameters in- 
between the separate fit results. We advocate separate description of strange hadrons 
and charmed hadrons. Improved data quality with large kinematic range (especially 
Q spectra at low px) "will provide better fit results for the multi-strange hadrons. A 
combination of spectra of J/ip and open charmed hadrons and/or their elliptic flow 
will enable us to break the correlation among /3, T and q. All these improvements 
await low-energy scan at RHIC, and future charm measurements at FAIR, RHIC and 

LHC [IHlEniEIlES]. 

This study also helps clarify the discrepancy on the claims about the early freeze-out 
of multi-strange hadron using a hydrodynamic evolution and a blast- wave model [5l[7j. 
In the Boltzmann-Gibbs blast-wave model (BGBW), the obtained /3 for bulk hadrons 
(7r,K and p) is about 0.6c, much larger than what we obtained from the TBW model 
{(3 ~ 0.5c). This may be due to the intrinsic assumption of BGBW where the local 
thermalization with a Boltzmann distribution increases the radial flow to absorb the 
non-equilibrium effect. The large /3 apparently increases the difference to that of multi- 
strange hadrons. On the other hand, it is not clear if freeze-out temperature is the only 
variable in the final hydrodynamic evolution, where the comparison of Q spectrum to 
the hydrodynamics was made [7] . Usually, it is argued that a hybrid model with hadron 
cascade at the hadronic phase can describe the experimental data better [53] . A least-x^ 
fit to the spectra using different freeze-out temperature settings in the hydrodynamic 
model will be valuable at achieving a quantitatively comparison among the different 
approaches. Our fit results in figure H] show strange hadrons freeze-out early but their 
radial flow is the same as that of light hadrons. This offers an alternative and novel 
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picture of hadron freeze-out: 

(i) At partonic phase, partons approach thermal equihbrium, resulting in small (g — 1) 
value for strange and light hadrons. However, ^/ip are not thermalized. 

(ii) After hadronization, multi-strange hadrons decouple from the system and their 
measured radial flow reflects that at this stage. 

(iii) At hadronic phase, hadron scattering doesn't produce a collective radial flow. It 
also is not sufficient to maintain the system in equilibrium. The consequence is 
that the copiously produced light hadrons are much further away from the thermal 
equilibrium at the end of the hadronic phase than at its beginning. This results in 
a large (g — 1) value for light hadrons without increasing the radial flow, this may 
be a consequence of the explosive nature of the fireball at early stage. 

(iv) This type of investigation is not possible with the BGBW since equilibrium is 
built in its assumption. However, TBW enables us to investigate the effect of 
non-equilibrium process. 

Recently, Tsallis statistics has been applied to extract chemical freeze-out 
information at RHIC and SPS [511 ES]- However, it shows that the resonance feed- 
down has a big impact on the minimum of the fit [55]. The chemical fit with 
Tsallis statistics doesn't produce a meaningful temperature parameter at the minimum 
value when the resonance feed-down is applied. It is tempting to conclude that 
resonance decays also have a big impact on the spectral shapes of stable hadrons at the 
kinetic freeze-out [7]. We argue that this is not the case. Models [56l [57], which allow 
the stable particles to continue the elastic interactions dominantly through resonances 
after chemical freeze-out, can explain the resonance measurements [T9l [20l [58] and the 
extracted different temperatures between chemical and blast- wave fits [H [5], [HI [IDl [33] . 
Although the grand total of each stable particle species including those from resonance 
decay at later time is defined to be conserved number after chemical freeze-out, for 
very short-lived resonances produced at chemical freeze-out, most (if not all) of them 
decay before the kinetic freeze-out. Meanwhile, their daughters from the decays interact 
with the other bulk hadrons in the medium and the correlation among the daughters 
as resonance from the chemical freeze-out has been erased. On the other hand, the 
stable particles and resonances after chemical freeze-out continue the elastic scattering 
and generate new resonances at a lower rate before kinetic freeze-out. The resonances 
available at the kinetic freeze-out are those regenerated from the bulk hadrons at the 
time close to the kinetic freeze-out. Therefore, the resonance yields and kinematics are 
completely determined by the kinematics of the daughter hadrons that generate those 
resonances, and not the other way around. As shown in [33], the effect of resonance 
decay has little impact on the final result. If any, it seems to support the idea that the 
resonances (i.e. p) are from regeneration with their yields and kinematics determined 
by the daughter kinematics. 
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4. Conclusions 

Tsallis Statistics in a blast-wave model (TBW) was used to investigate the non- 
Boltzmann distribution, radial flow and freeze-out temperature of particles and their 
dependence on particle species and beam energy in the relativistic heavy-ion collisions 
at SPS and RHIC. In addition to the light hadrons, ^/ip and the particles containing 
strangeness were examined separately to study their radial flow and freeze-out. The 
^ /ip m Pb+Pb collisions at SPS requires a temperature of 180 MeV and an average flow 
velocity of 0.06 c in the TBW flt, and therefore is consistent with early decoupling and 
obtains little radial flow. The J/^ spectra at RHIC are also inconsistent with the bulk 
flow proflle. However, there is a strong correlation among the obtained parameters from 
the TBW flt to the J/^ spectrum along. To improve the uncertainty and break the 
parameter correlations, higher statistics, elliptic flow and/or spectra of other charmed 
hadrons are needed to further determine the detailed J/?/' flow pattern at RHIC. We 
found that the strange hadrons approach equilibrium quickly from peripheral to central 
A+A collisions. The temperature from the TBW flt to strange hadrons is higher than 
that for light hadrons while their radial flows are similar. These results provide a novel 
picture of freeze-outs: a thermahzed system is produced at partonic phase; the hadronic 
scattering at later stage is not enough to maintain the system in equilibrium and does not 
increase the radial flow of the copiously produced light hadrons. The existence of such 
picture in reality seems to be natural. However, the investigation of such off-equilibrium 
effect has not been possible with the BGBW model or ideal hydrodynamics. 
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